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In a recent Letter [l| , Yepez et al. performed numerical 
simulations of the Gross-Pitaevskii equation (GPE) us- 
ing a novel unitary quantum algorithm with very high 
resolution. They claim to have found new power-law 
scalings for the incompressible kinetic energy spectrum: 
"...(the) solution clearly exhibits three power law regions 
for E^ orap (k): for small k the Kolmogorov fc~3 spec- 
trum while for high k a Kevlin wave spectrum of A: -3 ..." . 

In this comment we point out that the high wavenum- 
ber k~ 3 power-law observed by Yepez et al. is an arti- 
fact stemming from the definition of the kinetic energy 
spectra and is thus not directly related to a Kelvin wave 
cascade. Furthermore, we clarify a confusion about the 
wavenumber intervals on which Kolmogorov and Kelvin 
wave cascades are expected to take place. Finally, we 
point out that the incompressible kinetic energy spec- 
trum of the initial data chosen by Yepez et al. scales like 
at small wavenumbers, perhaps not the best choice 
to obtain a clean Kolmogorov regime. 

The dynamics of a superflow is described by the GPE 


(i) 


where the complex field ip is related by Madelung's trans- 
formation tp = y/pexp (i jt^ j to the density p and ve- 
locity v = \7(f> of the superfluid. In these formulae, £ is 
the coherence length and c is the velocity of sound (for a 
fluid of unit mean density) . The superflow is irrotational, 
except on the nodal lines ip = which are the superfluid 
vortices. 

The GPE dynamics Eq. (fT]) conserves the energy that 
can be written as the sum (the space-integral) of three 
parts: the kinetic energy = l/2( v /pu ) ) 2 , the inter- 
nal energy £ int = (c 2 /2)(p— l) 2 and the quantum en- 
ergy £ q = c 2 ^ 2 (dj^/p) 2 . Using Parseval's theorem, one 
can define corresponding energy spectra: e.g. the ki- 
netic energy spectrum Eki n {k) as the angle-average of 
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The 3D angle-averaged spectrum of a smooth isolated 
vortex line is known to be proportional to that of the 
2D axisymmetric vortex, an exact solution of Eq. (TTJ) 
given by ip vort (r) = y p(r) exp(±i</?) in polar coordinates 
(r, tp). The corresponding velocity field v(r) = \[2c^jr is 
azimuthal and the density profile, of characteristic spatial 

+ and \J p(r) = 


1 + 0(r~ 2 ) for r — > oo. Thus yfpvj has a small r singular 
behavior of the type r° and behaves as r -1 at large r. 
In general for a function scaling as g(r) ~ r s the (2D) 
Fourier transform is g(k) ~ k~ s ~ 2 and the associated 
spectrum scales as k~ 2s ~ 3 . Thus Eki n (k) scales as fc~ 3 
for k > k$ ~ £ _1 and as fc" 1 for k < 

Following the above discussion, the k~ 3 power-law ob- 
served in [l( is an artifact stemming from the definition 
of the kinetic energy spectra and is not directly related 
to a Kelvin wave cascade. 

Another very important scale, not discussed in the Let- 
ter p]], is the scale I of the mean intervortex distance. 
The energy cascade is expected to end at kg ~ £ _1 
and the Kelvin wave cascade to begin there, after an 
eventual bottleneck [4j. Note that £j ^> I ^ £, where 
li is the energy containing scale. We thus believe that 
nothing particularly interesting is taking place between 
k^ and the maximum wavenumber A: max of the simulation 
and that there is a confusion in [l[ between kg and k^. 

Furthermore, the initial data used in (see the sup- 
plementary material) is a 3D set of 12 straight vortex 
lines, with intervortex distance £ of the order of the box 
size. The A; -1 scaling of the initial data Eki n {k) thus ex- 
tends down to small wavenumbers kj ~ ij 1 . This behav- 
ior of the initial data is in contrast to the Taylor-Green 
initial data used in 3], where destructive interferences 
deplete the value of Eki n (k) in the interval ki < k < kp. 
The initial data chosen in pj is thus perhaps not the best 
choice to obtain a clean Kolmogorov regime. This might 
explain the high level of compressible kinetic energy in 
the —5/3 scaling zone that is apparent in the movie of 
the supplementary material. 

We thank M. Abid and C. Nore for useful discussions. 


extent £, verifies W p(r) ~ r as 
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